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ABSTRACT 

The formation of brown dwarfs via encounters between proto-stars has been confirmed with high- 
resolution numerical simulations with a restricted treatment of the thermal conditions. The new results 
indicate that young brown dwarfs (BDs) formed this way are disk-like and often reside in multiple 
systems. The newly- formed proto-BDs disks are up to 18 AU in size and spin rapidly making small- 
scale bipolar outflows, fragmentation and the possible formation of planetary companions likely as 
have recently been observed for BDs. The object masses range from 2 to 73 Jupiter masses, distributed 
in a manner consistent with the observed sub-stellar initial mass function. The simulations usually 
form multiple BDs on eccentric orbits about a star. One such system was hierarchical, a BD binary 
in orbit around a star, which may explain recently observed hierarchical systems. One third of the 
BDs were unbound after a few thousand years and interactions among orbiting BDs may eject more 
or add to the number of binaries. Improvements over prior work include resolution down to a Jupiter 
mass, self-consistent models of the vertical structure of the initial disks and careful attention to avoid 
artificial fragmentation. 

Subject headings: accretion, accretion disks — hydrodynamics — instabilities — stars: formation — 
stars: low-mass, brown dwarfs — stars: luminosity function, mass function 



1. INTRODUCTION 

Brown dwarfs (B Ds) are abundant, both in the field 
and in star clusters (|Chabrieill2?Kll iMartm et all 1200 It 
lLuhmanll2?KK? ). The fiat low mass end of observed stel- 
lar initial ma ss function (IMF) (e.^ ., Hamblv et a l. 1999; 
IMartm et al.| EM |Barrado y Navascues et al.| 
suggests that brown dwarfs are as numerous as hydrogen- 
burning stars. Observations indicate that a large per- 
centage of young BDs are surrounded by circumstellar 
disks (Muench et al. 2001). Mid-infrared, submillime- 
tre and millimetre measurements have characterized the 
spectral energy distributions (SEDs) of the BD disks 
(|Pascucci et al.l 12003). revealing the accretion rate and 
ongoing evolution of the young BD disks. The decreas- 
ing trend of IR excess in SEDs with age suggests that 
BD disks evolve from fiared to fiat geometry, similar to 
the evolution of disks surrounding large stars (|Mohantvl 
[2005). A T-Tauri phase has been det ected via broad 
Hey lines in many of these young disks (|Muzerolle et al.l 
|2003; Javawardhana et al. 2003)), some times associated 
with bipolar outfiows ( Whe lan et aTl 120051) . Giant 
planets have also be en imaged orbiting brown dwarfs 
dChauvin et aUHQQl. 

While the observational results admit the possibility 
that brown dwarfs share a common origin with stars (i.e. 
from the gravitational collapse of molecular cloud cores), 
the theory of the BD formation is far from well estab- 
lished. BD masses (0.013 — 0.075 Mq) are smaller than 
the typical Jeans mass in the parent clouds (~ 1 Mq) 
by more than an order of magnitude which seems to 
work against direct collapse of clumps as a formation 
mechanism. Approaches to solving this problem include 
lower ing the Jeans M ass via regions of low tempera- 
ture (|ElmegreenllT99l or via density enhancements in 
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I I200I ("turbulent 
fragmentation"). After low mass objects form ongoing 
mass accretion determines the final IMF. It has been pro- 
posed that all objects initially form with low masses and 
that observed BDs are those objects ejected from the gas 
cloud reservoir before they grow to stellar mass ("com- 
petitive accretion") (Bate et al. 2003; Reipurth & Clarke] 
I2OOID . 

Here we investigate a complementary mechanism of 
brown dwarf formation: instead of producing sub-stellar 
object in large numbers directly, we start by using the 
fact that initial cloud core-collapse stage forms numer- 
ous stellar mass objects with extended protostellar disks. 
Such early disks are very fiared and for this study they 
are assumed to be stable against spontaneous fragmen- 
tation. However, encounters between extended disks are 
very violent and can trigger gravitational instabilities 
and the formation of substellar objects. Since the gas 
density in a disk is much higher than an early-stage core, 
the local Jeans mass is comparable to a substellar mass, 
making it possible for gravitational collapse to produce 
objects that range down to planetary masses. 

The encount er pi cture has been studie d pr eviously by 
IWatkins et ^ ([T998a b ) and Lin e tHI (|T998) with low 
resolution simulations using proto-stars with fiat disks 
(two dimensional) or disks with highly-simplified vertical 
structure, respectively. The formation of substellar-mass 
objects was reported by both groups but it was not clear 
if the fragmentation was real or numerical or what the 
detailed properties of the objects might be. It is also 
worth noting that encounters occur in gaseous simula- 
tions of star cluster formation (e.g. IB ate et all ()2003l) ). 
However, modelling many hundreds of solar masses at 
once forces a compromise on the the mass resolution in 
individual proto-stellar disks and approximations (such 
as unusually low mass clusters or very large sink radii. 
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Fig. 1. — An encounter between two prograde rotating proto- 
stellar disks. The disk moving from top-right to bottom-left was 
tilted 7r/4 relative to the encounter plane. Panels 1-5 correspond 
to times 10000, 12000, 13000, 14000, and 15000 years, respectively. 
Panel 6 shows the fragmented disk at bottom-left of panel 5 with 
different zoom and density scale. Shock compressed gas (a), tidal 
structures (b) and strong disk instabilities (c) are indicated. The 
color bars indicate the density scale in logarithmic space in units 
of Mq/AU^. Note that the black regions in panels 1-5 are filled 
with gas, but have density < lO-^^ Mq/AU^. At 15000 years, BD 
Clumps are present in the central tidal feature and the inner disk 
(c) in the last panel. 

discussed in the next section) have been necessary. 

In this letter we report the results of new, high- 
resolution pairwise encounter simulations with realistic 
protostellar disk models and a range of initial configura- 
tions. We discuss the conditions under which fragmen- 
tation can occur and, in the cases where objects form, 
we examine their shapes, rotation and orbital properties 
and discuss their observational implications. 

2. SIMULATIONS 

Th e parallel TreeSPH code "Gasohne" (Wadslev et alj 
I2004D was used to simulate the encounters of protostellar 
disks. The initial condition s employed flared disk models 
(|Chiang Goldreichlll997D consistent with the observa- 
tions of early stellar objects. Each disk extended from 50 
to ^ 1500 AU and had mass ~ 0.6 M©, comparable to 
the mass of the central protostar (~ 0.5 M©). The radial 
surface density and temperature profiles were chosen to 
be similar to the intermediate quasi-static results of core- 
collapse simulations by Yorke & Bodenheimer (1999). 
The dependence is fiat inside 100 AU and then varies 
as r~^/^ until it falls off rapidly beyond 500-700 AU. 

The gas temperature was fixed according 

to the dista nce to the neares t protosta r as 

T(r) (X r~V^ (iMaver et al.H 2004: Pi ckett et all 1200,^ 
iChiang &: GoldreichI ll997Yand smoothly transitioned 
to the ambient temperature of ~ 40 K at r = 500 
AU. This temperatur e profi le is similar to that of 
lYorke Bodenheimer^ ("1999) which self-consistently 
included heating due to stellar irradiation, viscosity 
and shock waves which are beyond current simulation 
capabilities for non-axisymmetric simul ations. Based on 
Yorke's surface density and using the iD'Alessio et al. I 
(|2001 ) opacities the disks are optically thin outside 
clumps and cool efficiently. The product of the angular 
rotation rate and the cooling time ft tcooi = 0.3-0.7 
which is well within the required range for fragmen- 
tation for a molecular hydrogen disk of Q tcooi < 12 



iRice et al. I ElM UnhnsoTi ORmmie 1 12001. Given 
the importance of heating and the extremely short 
cooling time {tcooi ^ 30 years at 100 AU) it was felt that 
a fixed temperature distribution was most appropriate 
for the protostellar disks. Better treatment of the 
thermal conditions will be required for interior evolution 
of the brown dwarf clumps being pursued in future 
work. 

Vertical structure is one of the key factors affecting disk 
stability, lowering the critical Toomre Q valu e to values 
as low as ~ 0.6 for extremely ff ared disks (iKim et al. 
2002). Fo llowing arguments by lYorke Bodenheimer 
(1999) and lMaver et al.l(l2004D a nd experimentation with 
different O values (|Shen et al.l . in preparation), we se- 
lected an initial Qmin = 1.6 that is comfortably in the 
stable region for the simulations presented here. The 
disk vertical structure was constructed self-consistently 
assuming the gas to be in hydrostatic equilibrium with 
both its own global self-gravity and the star gravity. The 
scale heights in the midplane are quite large (~ 6 AU at 
the densest point at ~ 80 — 100 AU) and well resolved 
(local SPH smoothing 1.9 AU) and decrease above the 
mid-plane. The velocities were set to maintain circular 
orbits. 

Each disk was modeled with 200,000 SPH particles 
with masses of 3 x 10~^ (1 Earth Mass). An un- 
resol ved Jeans mass can re sult in artificial fragmenta- 
tion (|Bate BurkertI Il997t) . This resolution ensured 
that both the local Jeans mass and the vertical structure 
were resolved until well after any initial fragmentation 
that forms clumps. The disks were tested to be initially 
stable against spontaneous fragmentation. The artificial 
shock-capturing viscosity had no dynamical effect (mea- 
sured viscous timescales of 10^~^^ years). The encounter 
timescale (~ 10^ years) is too short for physical viscous 
effects to play an important role. 

The proto stars the stars were treated as sink particles 
((Bate et al.i2003 ) with radii of 1 AU but no new sink par- 
ticles were added during the runs. The softening length 
of gravitational force was chosen to be 0.2 AU, one tenth 
of the minimum initial spacing. Thus densities up to 
1000 times larger than the highest densities in the initial 
disk were well treated numerically. The softening was ex- 
pected to inhibit extremely low-mass fragmentation near 
the resolution limit. 

Pairs of disks were placed on initially unbound orbital 
trajectories with a projected impact parameter of 1000 
AU (prior to gravitational focusing). Various disk spin 
orientations (with respect to the angular momentum of 
the orbit) were run, covering prograde cases (spin-orbit 
angle less than 7r/2), retrograde cases (spin-orbit angle 
larger than 7r/2) and combinations of the two. 

3. RESULTS AND DISCUSSION 

The simulation results support the main conclusion of 
the previous studies: the formation of BDs can be trig- 
gered by protostellar disk encounters. However, the con- 
ditions for fragmentation are more restrictive than previ- 
ously reported, such that previously employed initial con- 
ditions (e.g. Qmin = 2.4) are unable to fragment when 
fully resolved. Vertically resolved disks closer to the criti- 
cal value of stability criterion are required {Qmin < 2.1). 
Convergence studies also indicate that artificial fragmen- 
tation is likely with resolution significantly less than that 
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Fig. 2. — Typical shape of newly formed Brown Dwarf mass 
objects. The specific object shown was produced via shock layer 
fragmentation in a coplanar disk encounter between disks rotating 
in a retrograde sense. The encounter occurred in the x-y plane. 
left: the object projected in the x-y plane; right: the same object 
projected in x-z plane. Both boxes are 3 AU x 3 AU across. 



employed for our production runs ( She n et all in prepa- 
ration). The primary new results, reported here, are the 
diskiness of the young BD, the extension of the mass dis- 
tribution down to planetary masses and the likelihood of 
multiple or even hierachical systems. 

A typical prograde encounter is depicted in Fig ^ 
Because of the tidal interaction, gas in the outer part 
{R > 300 — 500 AU) disperses. Fragmentation usually 
occurred after the disks passed periastron. In the case de- 
picted here, a 20 Mj^piter clump formed in a distant tidal 
arm structure at t = 14,000 years and two objects (with 
mass 9 and 70 Mj^piter) formed via gravitational insta- 
bility in an inner disk at t = 15,000 years. Three basic 
mechanisms were found for forming sub-stellar mass ob- 
jects: fragmentation within shock layers, condensations 
within prograde tidal structures and fragmentation via 
instabilities in the inner disks after inflows. Shocks (Fig. 
G] label a) formed due to the highly supersonic velocities 
of the impacts, but fragmentation in shocks only hap- 
pened in coplanar encounters. Tidal tails (Fig. ^ label 
b) were induced by gravitational resonance in prograde 
disks. Disk instability (Fig. ^ label c) was a common 
feature during encounters, as mass moved inwards due to 
angular momentum transport via tidal tails for prograde 
disks or cancellation of rotation in retrograde cases. 

The self-gravitating objects produced in the simula- 
tions are in the mass range from 2 (600 particles) to 
73 Jupiter masses, covering the BD and giant planet 
regimes. All of the objects are highly flattened and disk- 
like (Fig. A T-Tauri-like phase, implying a disk, has 
been observed for planetary mass object (Moha ntv et ahl 
120051) . At the end of the simulations, the proto-BD disks 
had radii 0.3-18 AU, with the lower limit quite likely re- 
lated to the gravitational softening (0.2 AU) and height- 
to-radius ratios of H/R ~ 0.1. The objects spin rapidly 
and most of the clumps formed via fragmentation of pro- 
tostellar disks have similar spin directions to the parent 
disk. But it is not necessary true for the clumps formed 
in shock layers or tidal tails. 

The simulations were halted several thousand years af- 
ter BD-mass objects form when the stars are well past 
peri-astron. During further evolution, gas in the BD 
disks with lower specific angular momentum should con- 
dense to form a central core, and the remaining disk will 
accrete onto the core on a viscous time scale. In order 
to for these fast-rotating disks lose their angular momen- 
tum, it is natural to expect observable bipolar outflows. 
These were confirmed observationally via the recent de- 
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Fig. 3. — A preliminary mass distribution of 32 objects resulting 
from 14 simulations. All proto-stellar disks in these simulations had 
similar vertical structure and were all stable against spontaneous 
fragmentation. The encounter velocities, the relative angle between 
two disks and the angle between the disk and the encounter plane 
differed from case to case. 

tection of several forbidden lines {e.g. [Oi]&l[Sii]) and a 
P Cygni-like dip in the iJ^v profile in the sp ectrum of a 
young BD p Oph 102 (|Whelan et al.lj200.^. The lower 
accretion rate expected in a BD disk (|Natta et al.ll200l 
makes it a scaled-down version of a T Tauri star as shown 
by the observed offset from the continuu m of the forbid - 
den emission lines in a BD-outflow ( Wh elan et al.ll2005D . 
Another way to distribute the angular momentum of the 
proto-BD disk is to form planets or further fragment. 
Thus BD companions can have planetary masses, form- 
ing a BD-planet system, or comparable masses to the pri- 
mary BD to form a binary or multiple BD system. The 
giant planet companion (2M1207 b) to a nearby young 
BD found by the new VLT/NACO imaging observations 
(|Chauvin et al.ll2(K)5l) is a possible example of the for- 
mer case. The large inferred inner holes in the spectra of 
some BD disks have be en conside red as signs of ongoing 
planet formation (Mohantvll2005 lD. 

Formation of binary BD pair occurred in one of our 
simulations via the secondary fragmentation of a large 
(6 AU) clump ~ 140 AU from the star. Though with 
5000 particles the clump is larger than the the en- 
tire protostellar disk in previous work, the Jeans Mass 
is only marginally resolved and this result should be 
confirmed with more resolution and additional physics. 
We note that fragmentation is not unexpected from 
a physical perspective given the high rate of spin 
iMatsumoto Hanawa I (|2003D . The two components (19 
and 13 Jupiter masses) orbited each other with separa- 
tion ~ 3 AU and orbital period ~ 30 years. Both retained 
a highly- flattened shape. It is likely that the system will 
remain as an hierarchical star-BD multiple. This case is 
observationally inferred for GL 569 where the primary is 
a M2.5V star (GL 569A) and a rec ently observationally- 
confirmed BD triple (GL 569B) (|Gorlova et all l20M 
iSimon e~ l. 2006). The age of the quadruple is about 
100-125 Myr so no disks are observed. However, in some 
of much younger BD multiples, accretion signatur es are 
indicated via a V-band excess (|Kraus et al.ll200d) . sug- 
gesting accretion disks around the components. 

The rapid gro wth assumption of the competitive accre- 
tion scenario (Reipurth & Clarke '2001"; Bate et al. 20031) 
makes ejection a requirement in order to retain BDs. In 
our simulations, about one third of the clumps (12 out of 
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32) were unbound and left the parent disks with ID RMS 
velocity dispersion v ~ 2.6 km/s, higher than the veloci- 
ties of the protostars. Among the bound BDs, three out 
twenty have large semi-major axes (a > 300 AU). As gas 
dispersion is severe in the outer disk region due to the 
interaction, these objects will lose access to a gas reser- 
voir and become free-floating BDs or wid e-separation 
BD companions as seen in observations fe.g.. lGizis et all 
1^01). For objects close to the star (a < 100 AU), ejec- 
tion may occur due to interactions with other orbiters. 
The need to invoke ejection to shut off rapid growth is 
still under debate. For example, iKrumholz et al.l (|2005D 
suggest that accretion is very slow and the final mass of 
an object is determined by the initial collapse and thus 
BD-mass objects should remain as brown dwarfs. The 
competitive accretion model also needs to better estab- 
lish that the first objects to form in a cluster environment 
will have such low masses. 

Encounter-induced protostellar disk fragmentation 
provides an additional mechanism to produce brown 
dwarfs and it does not rule out BD formation via other 
methods such as turbulent fragmentation. Ongoing work 
involving more encounter cases will build up statistics to 
enable estimates of which of these models is more likely 
to produce observed sub-stellar populations based on the 
production rates, BD disk properties, the frequency of 
multiples and the impact on proto-stellar evolution and 
accretion histories. 

A preliminary mass distribution of the simulated young 
BDs is plotted in Fig. The clumps were found 
to be more abundant at the lower mass end, which is 



broadly c onsistent with observed substellar initial mass 
function (|Martm et al.ll200ll iBarrado v Navascues et al.l 
I2004D . The turnover at the lowest (planetary) mass 
(M < 10 Mjupiter) is related to our resolution limits. 
An apparent gap in the population from 35-60 Mj^piter 
is suggested in the figure, but the current statistics are 
insufficient to confirm it and additional runs will be re- 
quired. Clumps typically form more than 100 AU from 
a star but sometimes significantly closer where the phys- 
ical conditions are more complex. However, there is no 
apparent bias (according to radius of formation) in the 
mass distribution and it looks similar if clumps within 
100 AU are excluded. 

If encounter-induced BD formation is dominant, the 
BD numbers (relative to the number of stars) is ex- 
pected to differ depending on the star-forming envi- 
ronment. Ob servation indications of this have been 
report ed fe.g ^Levine ^t al.| HHHl iLuhman et aQ HHHl 
Briceh o et al.ll2 002) but the data are not conclusive. The 
way this ratio varies as a function of cluster properties 
is not straightforward. For example, an increase in the 
cluster density indicates more impacts between its mem- 
bers, while a higher velocity dispersion in denser cluster 
greatly redu ces the number of fragments per encounter 
(Shen et al., in preparation). To make a concrete pre- 
diction regarding BD production as a function of star 
cluster properties, the more comprehensive set of sim- 
ulations now underway will have to be combined with 
detailed estimates of encounter rates in young star clus- 
ters. 
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